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ABSTRACT 

We have performed 2D hydrodynamic simulations of a pulsed astrophysical jet 
propagating through a medium that is populated with spherical inhomogeneities, 
or "clumps," which are smaller than the jet radius. The clumps are seen to affect 
the jet in several ways, such as impeding jet propagation and deflecting the jet 
off-axis. While there has been some debate as to the prevalence of these types of 
condensations in the ISM or in molecular clouds, the exploration of this region of 
parameter space nonetheless both shows the potential for these clumps to disrupt 
astrophysical jets and yields results which recover aspects of recent observations 
of Herbig-Haro objects. We find that the propagation of the jet and the vorticity 
induced in the clump/ambient medium correlate well with a "dynamic filling 
function" fd across all the simulations. 

Subject headings: ISM: jets and outflows - ISM: Herbig-Haro objects - hydrody- 
namics - ISM: clouds 



1. Introduction 

The investigation of astrophysical flows, such as those of Herbig-Haro objects, AGN 
jets, and planetary nebulae has benefited in recent years from the advancement of comput- 
ing capability. Hydrodynamic and magnetohydrodynamic simulations of HH objects have 
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been performed in 2D (e.g. Raga et al. 2002; O'Sullivan & Lery 2002) and 3D (e.g. de 
Gouveia Dal Pino 1999) for a variety of parameters. The effect of heterogeneous media on 
astrophysical flows has also received attention in terms of simulation studies. Adiabatic and 
non-adiabatic simulations of astrophysically-relevant shocked clouds have been performed 
by Poludnenko, Frank, & Blackman (2002) and Fragile et al. (2004), respectively, as have 
simulations of extragalactic radio jets in "clumpy" media (Saxton et al. 2005). Planetary 
nebulae simulations have also been carried forward which include the effect of clumps (Stef- 
fen & Lopez 2004). Typical motivations in these simulations are both the recovery of specific 
morphology and the capacity to characterize global characteristics, such as detailed velocity 
distributions. 

Observations of HH objects show them to be highly-collimated, several-hundred-AU 
radius structures with a primary bowshock 0.3 pc to several pc from their source with 
subsequent bowshocks or "knot" features trailing behind (see e.g. Hartigan et al. (2005) or 
the review by Reipurth & Bally (2001)). Detailed images of the bowshock reveal complex, 
dynamic structure. By comparing observations made 3 years apart, Bally et al. (2002) 
showed that the bowshock of HH 1 displays a wide range of motion in addition to its bulk 
flow. Using the two datasets, they were able to derive velocity vectors for the bowshock; see 
Fig. 1. The observations of multiple clumps in the flow made them particularly compelling 
though it was not clear if the clumps were fragments of the jet/bow shock disrupted via 
instabilities or environmental density perturbations enhanced by the passage of the shock. 

Of particular concern with HH objects is the relationship between the jet morphology 
and its environment. Jets which show large deflections have been argued to evolve from a 
precessing (Masciadri et al. 2002) or otherwise variable source (Cabrit & Raga 2000). Rosen 
& Smith (2004) performed pulsed jet simulations in which the density ratio between jet and 
ambient was varied. Earlier, Volker et al. (1999) performed a study in which pulsation, 
shear, precession, and opening angle were all varied. Further discussion of variable source 
studies may be found therein. 

While a variable source may indeed account for some "bending" of jets, there are other 
cases where the morphology demands an alternate explanation, such as HH 110. These 
explanations have typically included a significant interaction with an object in the jet's path 
(de Gouveia Dal Pino 1999; Raga et al. 2002) or with a strong ambient magnetic field (Dc 
Colle & Raga 2005). 

Thus the role of clumps influencing HH jets and HH objects is of interest for both 
observational and purely theoretical reasons. Jet interactions with heterogeneous media in 
other contexts such as planetary nebula are of interest as well for similar reasons. To date, 
simulations involving HH jets and clumps have considered inhomogeneities which are larger 
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Fig. 1. — 1994 epoch HST image of the HH 1 bowshock. The length of arrows corresponds 
to motion over 15 years. From Bally et al. (2002). 
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than the jet itself (see Raga & Bohm 1987; de Gouveia Dal Pino 1999). In this paper, we 
present results from 2D hydro dynamic simulations involving dense "clumps" which are 30 
AU in radius — a fraction of the jet's own 100 AU radius. 

The existence of dense condensations on the order of 10-100 AU has been the subject 
of some debate. In the ISM, Heiles (1997) asserted that features which appeared to be 
these small clumps could in fact be larger structures — filaments or sheets — which are viewed 
edge-on. Deshpande (2000) suggested that the supposed observations of clumps of this size 
could in fact be accounted for by a power law extension of irregularities in opacities on larger 
scales in the ISM. However, Faison & Goss (2001) presented results using very long baseline 
interferometry observations of several extragalactic sources which they purported to show 
several distinct clumps less than 100 AU in extent. More recently, Hennebelle & Passot 
(2006) presented work arguing that the presence of Alfven waves may serve to enhance the 
formation of small-scale structure in the ISM. 

In star- forming regions, detection limits play an increased role. Vannier et al. (2001) 
observed H 2 from a region of the Orion molecular cloud with a resolution of ~0.15", cor- 
responding to physical length of ~50 AU. However, synthetic data used for calibration of 
Fourier analysis created spurious signals below ~100 AU, and so data below this level were 
not used. With regards to the prevalence of small-scale structure in star-forming regions, the 
situation seems even less clear. The survey of Hatchell et al. (2005) of the Perseus molecular 
cloud reported beam mass sensitivity at 12 K only down to O.4M — the structures here in 
question are considerably smaller. 

At the very least, the observation and explanation of media on this small scale has 
proven to be difficult. While the existence of clumps on this scale remains an open question, 
the investigation in this area of the parameter regime remains of interest. As we shall see 
the propagation of radiative jet in a clumpy media yields some new results which will be of 
general interest. 

2. Computational Method & Physical Model 

We here employ the 2D (slab-symmetric) hydrodynamic, radiatively-cooling AMR ca- 
pability of the AstroBEAR computational code 1 . The numerical code solves the conserved 
form of the equations of fluid dynamics with radiative cooling provided by the cooling curve 



1 Information about the BE ARCLAW / AstroB E AR code may be found online, at 
http: / /www. pas. rochester.edu/~bearclaw/ 
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of Dalgarno & McCray (1972) and a monotonic ideal gas equation of state. The AstroBEAR 
code has been tested on a variety of problems (see e.g. Poludnenko, Frank, & Blackman 
2002; Varniere et al. 2005; Cunningham, Frank, & Blackman 2006). 

Six runs altogether were undertaken, as discussed below. The resolution of the simula- 
tions were 312 by 156 on the coarsest level with three levels of refinement (for a maximum 
resolution of 2496 by 1248). This resulted in 52 and 16 cells per jet and clump radius, 
respectively, on the finest level. Running on 8-12 processors, each simulation took roughly 
three days to complete. All simulations were allowed to run until the jet propagated off the 
end of the domain, which occurred between t=162-554 years, depending on the simulation. 

We have used the observations of Bally et al. (2002) as motivation for our investigation. 
In our simulations, a jet 100 AU in radius with number density rij = 10 3 cm -3 and tem- 
perature Tj = 10 4 K propagates into an ambient region 4800 by 2400 AU in extent having 
number density n a = 100 cm~ 3 and temperature T a = 2, 000 K. The jet source is varied 
sinusoidally, 



where Vj = 150 km/s, v a = 75 km/s, and r = 30 years. This velocity modulation creates a 
series of internal pulsed bowshocks along the jet, as expected (Raga & Kofman 1992). 

We introduce dense circular "clumps" into the environment with radii of 30 AU and 
temperature set so as to maintain pressure balance with the ambient medium. The clouds 
are all the same size and have a hyperbolic tangent density profile, 



where the inclusion of the 0.1r c factor serves to smooth the outer 10% of the clump, keeping 
most of the clump area at a constant density n c . The number of clumps and their density 
are varied in the six runs — see Table 1. 

The positions of the clumps were seeded using random number generation with the 
same seed for each simulation. This allowed for continuity across simulations, as those 
clumps present for iV c = 60 exist also for N c = 120, and those at N c = 120 carry over 
to N c = 240. A selection algorithm was applied so that the minimum allowed separation 
between clumps would be 2r c . The clumps were chosen to fill the restricted domain (x=500- 
4800 AU, y=0-2400 AU) so that the jet head would have time to form before impacting the 
condensations. 



v(t) = Vj + v a sin 



(1) 
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Passive tracers are used to identify both clump and jet material in the simulation. This 
proves to be most useful in the analysis of vorticity and material mixing, as discussed in § 4. 

In figures comparing the six runs A-F, we present them in approximate order of increas- 
ing clump influence. The baseline case of a jet propagating in a smooth medium without any 
clumps is Run A. Defining the case of 120 clumps with a density ratio between clumps and 
the jet of n c /nj = 100 as the "standard" setup, we investigate the effect of clump density in 
two runs, with density ratios n c /rij = 1 and n c /rij = 10 denoted as Run B & C, respectively. 
The clump baseline run itself is Run E. Finally, we investigate the effect of varying filling 
factor by halving and doubling the number of clumps in the domain in Runs D (60 clumps) 
and F (240 clumps), respectively. 

We define the "spatial filling factor" as the ratio of the total area of clumps to that of 
the ambient environment, / = (N c Trr^/A a ), where r c and N c are clump radius and number, 
respectively, as given in Table 1, and A a is the ambient area. We further define the "dynamic 
filling factor" as the clump-to-ambient mass ratio, fd = (p c /Pa)f- 

3. Description of Simulations 

Figure 2 gives the results for the simulations at the same time, t = 180 years, with the 
images ordered top-to-bottom as in Table 1. 

1. The jet-only case (with the jet already propagated off the grid at this time) shows a 
smooth bowshock with some disruption due to Kelvin-Helmholtz instabilities (Blondin, 
Fryxell, & Konigl 1990). Internal working surfaces are coherent. 

2. The next image shows the jet with 120 clumps of density ratio n c /rij = 1. The 
basic shape of the bowshock is the same as in A), but clumps now "pockmark" the 
surface, causing slight localized disruptions. Internal working surfaces, and the jet's 



Run n c 


cm 3 


njrij 


N c 


/ 


fd 


A.. 












B.. 


10 3 


1 


120 


0.03 


0.295 


C. 


10 4 


10 


120 


0.03 


2.95 


D.. 


10 5 


100 


60 


0.015 


14.7 


E.. 


10 5 


100 


120 


0.03 


29.5 


F.. 


10 5 


100 


240 


0.06 


58.9 



Table 1: Parameters for the six simulations A-F. See text for details. 




Fig. 2. — Synthetic Schlieren image of Runs A-F all at the same time in the simulation. See 
text for details. 
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propagation, are essentially unchanged. The conclusion is that clumps at this density 
ratio are too light to noticeably affect the jet. 

3. The third image has an increased density ratio of n c /rij = 10. Here, the bowshock is 
more clearly disrupted. While the bowshock feature near (x = 2200 AU, y = 1900 
AU) is still present, other features of the jet's K-H rollovers have been obscured by the 
clumps. Reverse-facing bowshocks around clumps are present, as is the disruption of 
the first internal working surface. For the first time, the propagation of the jet has been 
retarded. The head of the jet is also wider than before, indicative of its disruption. 

4. Now the density ratio has increased to its maximum value of n c /rij = 100, but the 
number of clumps has been halved, from 120 to 60. While the overall shape of the 
bowshock is similar to the jet-only case of A), the clumps are seen to have great 
influence in their local neighborhoods. Prominent reverse-facing bowshocks are seen, 
as is the disruption of the internal working surfaces beyond a certain point (roughly 
x = 1800 AU at this time in the simulation). The width of the bowshock near the 
leading edge is again widened as compared to Run A. However, its propagation is not 
as slowed as in the previous image, suggesting a greater dependence on filling factor 
than density ratio — see the next section for a quantitative assessment of this effect. 

5. The density ratio remains the same, and the number of clumps is brought back to 
120. The effect of the clumps both on propagation and morphology is great. Internal 
working surfaces are disrupted to the degree that they are "piling up" near x = 2000 
AU at this point in the simulation. The clumps have also served to deflect some of 
the jet's momentum radially (upwards): the feature at (x = 3200 AU, y = 1700 AU), 
engendered by clumps in the range (x = 1800 — 2400 AU, y ~1100 AU), will continue 
to travel at an angle ~30° from horizontal. The width near the jet head has increased 
substantially, and the propagation has been reduced by ~38% from Run A. 

6. Finally, the number of clumps is doubled from 120 to 240. Here, the jet's propagation, 
disruption, and deflection are all vastly different from the jet-only case of Run A. A 
large part of the jet is being deflected upwards at an angle of ~25°, and the jet is 
dominated by small shocks around the clumps. The propagation of the jet has been 
reduced significantly, by roughly 58%. 

As can be seen in the bottom two panels of Fig. 2, the clumps can indeed have a significant 
impact on deflecting the jet, though it requires high density ratios and numbers. Whether 
or not this jet deflection is a transient phenomenon depends highly on the particulars of 
the clump distribution. Unless the distribution displays a high degree of coherent structure, 
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subsequent pulses of the jet will eventually bore a hole through the clumps, sweeping them 
to the side or entraining them. Eventually, the jet will reestablish flow primarily along the 
axis of the jet. This effect would probably change a great deal in 3D geometry, as discussed 
in § 5. 

To illustrate the evolution of the jet as it propagates through the inhomogeneous 
medium, we present Fig. 3. Here, Run E (n c /rij = 100, N c = 120) is shown at six times 
throughout the simulation, every ~90 years from t — to 440 years. As the jet pulse has a 
period of 30 years, the internal working surfaces appear to be at the same location throughout 
the images, while in truth they capture every third pulse. The panels, from top-to-bottom, 
show 

1. The initial clump-filled domain. (Note that a description of the clump-seeding algo- 
rithm is given in the preceding section.) 

2. The jet has already propagated halfway across the domain. Several clumps near the 
jet axis are being overrun by the jet's bowshock, though only a few reverse shocks have 
formed. 

3. The jet propagation has slowed, largely due to momentum being redirected at an angle 
of ~ 30° from horizontal. A consequence of the deflection is the stagnation of material 
on the axis, which will be swept away between this and the subsequent image. 

4. The head of the jet has largely reclaimed the leading position, though the effects 
of deflection both above- and below-axis remain evident in the width and angle of 
bowshocks on the wings. Note also the "swirling" motion near the base of the jet, due 
to the complex interaction between the jet launching and nearby clumps. 

5. The leading edge of the jet has propagated off the domain, and the jet's bowshock 
continues to expand radially. With few exceptions, no clumps which were originally 
within the jet's radius now remain. Additionally, the interior of the jet cocoon is almost 
entirely dominated by reverse-facing bowshocks. 

6. Finally, in the last frame, most clumps which still exist possess well-developed bow- 
shocks. Many of those near x = 1000 AU exhibit vertical "streaming" features whereas 
those downstream have broader shapes. All clumps within the jet radius have been 
swept clear, and as a result the internal working surfaces are beginning to recover 
coherent structure. 




Fig. 3. — Synthetic Schlieren image of Run E at six times throughout the simulation. See 
text for details. 
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4. Analysis 

The propagation of the jet through the domain may be investigated by tracking the 
leading edge of the jet throughout the simulation, as shown in Fig. 4. Here again the effect 
of the clumps on the jet is clear, as is the fact that the number of clumps has a far greater 
effect than the density thereof. The crossing times in years for Runs A-F are, respectively, 

t end = (162, 165, 202, 180, 277, 554) . (3) 

In the jet-only case, this corresponds to a bowshock speed of ~30 AU/year. This may be 
compared to the theoretical bowshock speed of Blondin, Fryxell, & Konigl (1990), 



v bs - V -j== (4) 

which for Vj = 150 km/s gives a speed of Vb s — 32 AU/year and a crossing time of 150 years, 
in fair agreement. 

The variability of the speed of the leading edge between simulations suggests two natural 
perspectives with which to discuss their evolution — namely, time and the position of the 
leading edge of the jet. Comparisons involving the former compares processes occurring at 
the same absolute time in the simulation, the latter serves to effectively describe them on 
the same scale of jet propagation through the medium. While both are legitimate, the latter 
perspective lends itself to a more illuminating investigation of the simulation suite in nearly 
all details. The investigation of vorticity benefits from both descriptions, as will be seen. 




Fig. 4. — The leading edge of the jet material as a function of time for each run. Symbols 
for each run are Run A: solid, Run B: dashed, Run C: dot-dashed, Run D: dotted, Run E: 
+, and Run F: o. 

Using a nonlinear least-squares fit, the time needed for the jet to reach the edge of the 
grid is relatively sensitive to the clump number {t en d ~ A" C L5 ), and is only weakly dependent 
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on their density (t end ~ n° 31 ). Furthermore, the dependence of grid-crossing time on the 
factor f d is approximately squared (t end ~ /J' 05 ); see Fig. 5 . 




Fig. 5. — A power law fit to the crossing time t end as a function of the dynamic filling factor 
f d for Runs B-F. The power law fit corresponds to an index of 2.05. 

With such a complex interaction between jet and clumps, the questions of mass loading 
and momentum entrainment come to the fore. Namely, is all the clump material swept to 
the sides by the jet's bowshock or is there significant bulk axial momentum imparted to the 
clumps? How much and in what manner does the jet and clump material mix? 

By introducing the clump filling factor f^, similar to / but for a given run R, and 
the jet filling factor fj* (analogously), we may investigate the evolution of their respective 
mixing and diffusion as the jet propagates through the grid. Figure 6a & b show and f^, 
respectively, versus the position of the jet leading edge. The mixing is seen to be coupled 
with the dynamic filling factor f d . By the time the jet reaches the end of the domain, the f c 
value for Runs B-D is only ~0.25, but for E and F is 0.47 and 0.77, respectively. This stands 
to reason as the slower propagation speeds in Runs E and F allow more time for the clumps 
behind the relatively stalled jet head to be disrupted into the environment. Similarly, Panel 
b demonstrates that the jet's lateral expansion is enhanced by increased f d . By the time 
the jet has reached the edge of the grid, the jet filling factor for Run F is more than twice 
that of Run A (// = 0.35 and ff = 0.84). The measure of f d thus provides a quantitative 
description of the disruption of the jet by the clumps: though its forward progress is slowed, 
ultimately the jet will increasingly fragment and fill the domain as it progresses through the 
medium. 

The fact that the jet and clump filling factors depend similarly on f d suggests that 
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Fig. 6. — The evolution of the clump filling factor f c and jet filling factor fj as functions of 
the position of the leading edge of the jet. The top panel is a, the bottom b. Line styles are 
as in Fig. 4. 
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Fig. 7. — Clump filling factor f c plotted as a function of jet filling factor fj for the five runs 
B-F. Line styles are as in Fig. 4. 
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there may be a correspondence between fa and the coupling between jet and clump material. 
Figure 7 demonstrates this: there exists a higher correlation — defined by the slope of the 
linear trends of the data — between f c and fj for higher As fd increases, the correlation 
number relating the two quantities rise from ~0.35 for Run B to ~0.80 for Runs D-F. This 
is to be expected for as jet material is deflected by the clumps, clump material is spread by 
the jet. The mixing of both materials are indeed dynamically coupled. 
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Fig. 8. — Binned and normalized speed of jet and clump material for Run E at the same 
time as in Fig. 2. Domain is binned into 12 sections along y, with velocity summed over the 
entire x-extent for each bin, then normalized over the collection. Jet material is solid line 
and clump material dashed. 

Momentum entrainment also appears to occur, as shown in Fig. 8. In this figure we 
divide the radial extent (i.e., along the y-axis) into twelve bins. Within each bin we sum 
over the entire axial extent of the grid (i.e., over all x) the axial component of the velocity, 
and normalize the results. Then, the bin with the highest velocity may be compared with 
bins which show less motion. Here we see that the further out one progresses radially, the 
less the axial momentum is transferred to the clumps. These data may be compared with 
Fig. 17 of Bally et al. (2002), which analogously gives the proper motion of features as a 
function of distance from the flow axis. The figures are in fair agreement. 

The possibilities of turbulent clump mixing is addressed by Fig. 9. The total magnitude 
of the vorticity of all material in the grid is dominated by processes associated with the jet 
flow (e.g., the K-H rollover). In lieu of investigating this total vorticity u>t, we investigate 
the evolution of the fractional clump vorticity uj c /ujt as determined by the clump tracer 
versus time and leading edge position in Fig. 9. The top panel, in conjunction with the 
crossing-times listed in Eq. 3, demonstrates that the fractional vorticity flattens (Runs D & 
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Fig. 9. — The fractional clump vorticity u) c /u T as a function of time and position of the 
leading edge of the jet for runs B-F. Line styles are as in Fig. 4. 




Fig. 10. — A power law fit to the fractional clump vorticity uo c /uot as a function of the clump 
filling factor f c for Run F. The power law fit corresponds to an index of -1.1. 
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E) or decreases (Runs B & C) once the jet propagates off the grid. Although not depicted on 
the plot, Run F eventually flattens to a value near that of Run E. Further, the rate of falloff 
varies inversely with f^, suggesting that the ability of an environment to sustain vorticity is 
directly correlated to the clump densities within that environment. 

In the bottom panel, two effects are seen. The first is that increasing clump density n c 
affects the scale of fractional vorticity in a regular manner: the increase in final u c / oot from 
Run B to C and Run C to E are the same, showing a consistent increase of 11% between 
each pair. Secondly, an increased clump number N c affects the concavity of the fractional 
vorticity's evolution. That is, while the final values are nearly the same for Runs D-F, the 
rate at which this value is attained increases markedly with number N c . 

Additionally, the dependence of fractional clump vorticity iv c / ujt on clump filling factor 
f c has a power law index of roughly -1.1. This relation is shown in Fig. 10 for Run F. 
However, it holds for all Runs B-F, indicating that the generation of turbulent motion and 
clump dispersal are well correlated. 

Finally, as a purely qualitative result we present Fig. 11, which is a snapshot of the 
standard run with an array of velocity vectors overplotted for regions where the clump tracer 
and velocity are nonzero. This may be compared with Fig. 1 and some basic similarities can 
be seen. Note that dispersal of velocity vectors near the head showing the change in direction 
of the jet. The redirection of the flow near clumps is also apparent. Additionally, vectors 
further from the jet axis are roughly half or less the magnitude of those near the axis as 
discussed previously. Of further note are the velocity vectors near the jet launching region 
which are directed perpendicularly or even backwards, indicative of the vortical motions 
present within that region. 

5. Discussion and Conclusion 

We have presented simulations of a pulsed astrophysical jet propagating into a hetero- 
geneous medium which possesses small-scale condensations. We find that the condensations 
can have an effect on the jet propagation which is strongly correlated to the filling factor of 
the distribution and less so on the density contrast between the jet and clumps. 

It is reasonable to anticipate that the slowing of the jet will allow increased time to mix 
clumps, resulting in increased filling factors for the clumps. If the clumps and the jet are 
dynamically coupled, this will result in an increased jet filling factor as well. Indeed, when 
considered against the position of the leading edge, both quantities are noticeably enhanced 
with increased fa- 
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Fig. 11. — A synthetic Schlieren image of Run E at the same time as in Fig. 2, now with 
velocity vectors overlaid. Details of note include large off-axis vectors near the leading 
edge and reverse-pointing vectors near the base of the jet due to vorticity. A superficial 
correspondence may be seen between this and Fig. 1. 



-18- 



Similarly, the fraction of the vorticity in the clumps has a clear dependence on the dy- 
namic filling factor when considered against the position of the jet's leading edge: increasing 
fd results both in a higher fraction of the vorticity attributed to the clumps and in the earlier 
realization of this higher fraction. It may therefore be said that fd is a robust indicator of 
the generation and sustainability of vorticity in the domain. 

Qualitatively our results can be considered in the broader context of astrophysical jet 
research. In particular, the morphologies presented here bear resemblance to those of Saxton 
et al. (2005), who performed 2D relativistic hydrodynamic simulations of overpressured, 
underdense AGN jets propagating into a heterogeneous medium. Unlike the present work, 
the heterogeneity was seeded via a Fourier algorithm intended to mimic power-law structure 
formation by turbulence, resulting in a variety of shapes and sizes of inhomogeneities. Despite 
this disparity in the scenarios, both simulations shared many morphological features, perhaps 
suggestive of scale-robust physics. See for example Figure 6 of that work, in which a jet is 
deflected off-axis in much the same manner as our standard run. If we cast those runs in 
terms of the dynamic filling factor fd as defined here, their simulations take on values 4.11- 
118.0, covering some of the same range as in the present study. While one might be tempted 
to assume the highest correlation between the two studies would occur for runs with the 
most-similar fd, the variable sizes of inhomogeneities in the previous work prevent such a 
strict correspondence to be drawn. In fact, the two simulations most similar are that study's 
Run Al and the present Run E, with values for f d of 14.6 and 29.5, respectively. 

Both of the results given in Figures 7 and 8 suggest that viable mass loading and 
momentum entrainment can occur. This is in contrast to the mechanism discussed by Raga 
et al. (2003). In that case, an object impinged upon the side of a steady propagating jet and 
was found to entrain only transiently. One possible explanation (as mentioned in that work) 
is that whereas the contribution from a single clump would have little long-term effect, the 
continued entrainment by successively over-run clumps serves to maintain the mass loading. 

Another result from our study is the realization that objects on size scales r c < rj 
can have a large influence on an astrophysical jets. Our simulations show clumps can of- 
fer a significant disruption mechanism in radiative (stellar) jets which has not before been 
considered. 

As noted however in the preceding section, many of the results presented here would 
doubtless be different in a 3D geometry. In particular, one could imagine that the impediment 
presented by the clumps would be diminished due to the additional degree of freedom. This 
suggests that the recovery of similar results would require either higher density ratios or 
larger filling factors. 
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Additionally, magnetic fields appear to have an effect on jet propagation (see e.g. Frank 
et al. 1998; Gardiner et al. 2000; Gardiner & Frank 2001; De Colle & Raga 2005), and the 
fact that magnetic confinement has been advanced as a possibility for maintaining small- 
scale structure (Diamond et al. 1989) suggests MHD simulations in this parameter regime 
may be useful. 

We now consider the observational evidence for structures on these scales beginning 
for reference with the ISM and then examining molecular clouds. As we shall see the issue 
remains contentious. The presence of structure on 10-100 AU scale in either the ISM or 
star-forming regions is not yet settled. 

Heiles (1997) demonstrated the difficulty of having dense (10 4 — 10 5 cm -3 ), small-scale 
(~15 AU) atomic structure be in pressure balance with the ISM due to the necessarily low 
temperatures involved. Of the three cooling mechanisms mentioned — the absence of dust 
grains, cooling via CO refrigeration, and an attenuated interstellar radiation field — only 
CO cooling was considered reasonable to entertain, though it was still unsatisfying. As an 
alternative, he put forth the idea that apparent small-scale features are in fact edge-viewed 
larger structures such as filaments or sheets. Heiles demonstrated that the requirement on 
pressure balancing is ameliorated somewhat by this concept, and discussed the observational 
consequences. (However, this conclusion does not directly negatively affect our results, as a 
2D slab-symmetric simulation mimics exactly these types of structures.) 

Subsequently, Deshpande (2000) presented analytic results purporting that small-scale 
structure in the ISM does not need to exist. Instead, irregularities in large-scale structure 
can mimic small-scale structure. This conclusion was reached by extending the power-law 
spectrum of structure on scales ~0.1 pc down to 10 AU. 

However, Faison & Goss (2001) presented new observations which seemed to contest 
this assertion. Using very long baseline interferometry, Faison and Goss were able to resolve 
absorbing features of order 20 milliarcseconds. Using 21cm absorption of galactic H I toward 
several bright extragalactic continuum sources, they were able to resolve several distinct ob- 
jects with inferred densities ~10 4 cm -3 and scales ~10-100 AU. Of the four sources imaged, 
three returned results which appeared in line with Deshpande's treatment (a variation in 
optical depth of about 0.1). One, however, showed a change in opacity of 0.6-0.8, quite in 
excess of the prediction. Moreover, Deshpande's analysis suggested that only gradients in 
optical depth should be observed since the variations were coming from larger-scale struc- 
ture. However, in each case presented by Faison and Goss, distinct objects were apparent, 
in direct contradiction to that prediction. 

Relatedly, Lauroesch, Meyer, & Blades (2000) have investigated small-scale structure 
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via Na I absorption toward a binary star system. While small-scale structure did seem 
present, the densities were significantly lower than those typically claimed (only ~20-200 
cm -3 ). 

For star-forming regions, the presence and prevalence of structure on this scale is less 
clear. For the densities and temperatures used here, emission would be expected either 
in molecular lines (for the low-contrast clumps) or lines more reminiscent of cloud cores 
(for the heavier clumps) (Vannier et al. 2001). In an investigation of structure sizes in the 
Orion molecular cloud, Vannier et al. (2001) found a preferential size scale between 3 • 10~ 3 
and 4 • 10~ 3 pc (600 and 800 AU) via area-perimeter, Fourier, and brightness distribution 
analysis of H2. However, their study was biased against small-scale structure due to a limiting 
resolution of 0.15", corresponding to a physical length of 50 AU. In particular, the model 
against which the Fourier data was vetted showed spurious small-scale signals due to noise 
below ~5-10~ 4 pc (100 AU) with a maximum at 3.5 • 10~ 4 pc (75 AU). As this peak was 
also featured in the real data, these and smaller scales were not included. Attempting to 
directly measure and count structures was once again impacted by the resolution of 0.15". 
These results suggest that structure below these scales in Orion is at least ambiguous, if not 
wholly endorsed. 

As to prevalence, Hatchell et al. (2005) imaged the complete Perseus molecular cloud at 
850 and 450 /im in order to map locations and conditions for star formation. However, their 
study surveyed cloud cores only down to a mass of O.4M . Assuming spherical symmetry, a 
single clump in the present work is substantially smaller than this limit. The Galgarone, Hily- 
Blant, & Levrier (2004) review noted that, considering the positive correlation between gas 
mass and size scale, and assuming that this trend holds for all scales, then it is reasonable to 
assume that small-scale structure cannot comprise a notable fraction of the molecular cloud 
mass. This does not however preclude their presence in regions close to cloud cores where 
stars have recently formed and where jet propagation is initiated. 

It should be noted that pressure balance between the clumps and ambient medium 
forces the clump temperatures to the range 2-200 K. Similar overall parameters have been 
simulated in prior studies; in many cases, pressure balance appears to have been sacrificed to 
avoid egregiously cold clumps (e.g. Raga et al. 2003; Melioli, de Gouveia Dal Pino, & Raga 
2005). While the simulations are most likely robust to the exact details of the temperatures, 
this aspect of the parameter regime deserves further careful consideration. 

Our conclusion is that, though the presence of small clumps remains uncertain, our 
results address significant dynamics which have heretofore been unexplored. Jets clearly 
show clumpy morphology and we have explored one means by which these morphologies 
can be achieved and the consequences of propagation in heterogeneous media. We can also 
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suggest a new direction in which to turn in consideration of clumpy jets. While clumps may 
form from instabilities associated with cooling at the jet head, its is also possible the jet 
beams are themselves clumpy. Future studies will explore distinct clumps along the axis in 
the beam which interact internally with the jet. It is reasonable to suppose that a clump so 
launched could overtake a bowshock within the jet and disrupt it. The exact details of this 
interaction should be a matter of some interest. 
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